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Abstract: We present a next-to- leading order (NLO) treatment of the production of a 
new charged heavy vector boson, generically cahed W, at hadron cohiders via the Dreh- 
Yan process. We fuhy consider the interference effects with the Standard Model W boson 
and allow for arbitrary chiral couplings to quarks and leptons. We present results at both 
leading order (LO) and NLO in QCD using the MC@NLO/Herwig++ and POWHEG methods. 
We derive theoretical observation curves on the mass-width plane for both the LO and 
NLO cases at different collider luminosities. The event generator used, Wpnlo, is fully 
customisable and publicly available. 
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1. Introduction 

There exists a proliferation of theories which contain new heavy, neutral or charged, gauge 
bosons referred to as Z' and W respectively. Both the Z' and W have been studied 
extensively and recent reviews can be found in [1] and [2]. The present study focuses on 
W bosons. The new charged vector bosons may or may not have similar properties to 
the SM bosons, depending on the theory that predicts them [1]. In particular they may 
have right-handed instead of left-handed couplings, may couple to new fermions, or may 
even be fermiophobic. Popular models which predict new charged vector bosons are based 
on extensions of the electroweak gauge group, SU{2) x ?7(1), for example to the gauge 
group SU {2)i X SU {2)2 x U{1) [3], or groups that contain the electroweak symmetry, such 
as SU{2>) X U{1) or SU{4) x U{1) [4]. Several models with extra dimensions contain W' 
bosons as Kaluza-Klein excitations in the bulk. Examples of these models include the 
Randall-Sundrum model with bulk gauge fields [5] and Universal Extra Dimensions [6,7]. 
Theories which break the electroweak sector dynamically may also contain the W as a 
composite particle [8,9]. 
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Current Monte Carlo simulations of Drell-Yan W' production at hadron colliders rely 
mainly on leading order QCD matrix elements and parton showers. There currently exists 
no treatment of next-to-leading (NLO) QCD effects which simultaneously includes the full 
interference effects for the W. In the present paper, we present the results of the event 
generator package Wpnlo [10] which improves the treatment of Drell-Yan production of 
heavy charged gauge bosons. We consider the interference effects with the Standard Model 
W, which have been shown to provide valuable information [11], but have not yet been 
considered in experimental searches. We use the 'Monte Carlo at Next-to-leading Order' 
method [12,13] for the Herwig+-|- general purpose event generator [14] (MC@NLO/Herwig++) 
and the 'Positive Weighted Hardest Emission Generation' method (POWHEG) [15, 16] to 
match the NLO QCD calculation to the parton shower, producing fully exclusive events. 
Note that a similar implementation of the Z' exists for the NLO MCSNLO event generator, 
which matches the complete NLO matrix elements with the parton shower and cluster 
hadronization model of the Fortran HERWIG event generator [17]. 

The paper is organised as follows: In Section ^ the leading order reference model is 
presented in detail, including the relevant assumptions, and the W — W' interference effects 
are studied. In Section | the MC@NLO/Herwig++ and POWHEG NLO implementations of the 
W' reference model are discussed. Section § presents a sample of distributions obtained 
using the Wpnlo event generator package. In Section ^ we present a theoretical discussion of 
the extraction of mass- width observation limits for the W at LO and NLO. We present our 
conclusions in Section ^. The appendix contains the description of a theoretical analysis 
for discriminating between models. 

2. W at leading order 
2.1 The reference model 




Figure 1: Feynman diagram for pp — > W/W — > li^X. 

The W' reference model is based on the one which originally appeared in [18]. In 
the model described therein, sometimes referred to as the Sequential Standard Model, the 
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W' couplings to fermions are directly transcribed from the SM W, i.e. it is a heavy copy 
of the SM W. In the present treatment we allow both right- and left-handed couplings, 
corresponding to (1 ±75) respectively, as well as arbitrary mixtures of the two. In the case 
of right-handed couplings, we assume that the right-handed neutrinos are light compared 
to the W, but not light enough for the Z boson to decay into them. The W and W 
couplings to fermions are given by (for i = W, W ) : 



V V2 



2 \ 1/2 

J - h^J,)f'Wt + h.c. (2.1) 



where Gp is the Fermi coupling constant, M\y is the SM W mass, C^''^ are the coupling 
strengths of boson i to leptons and quarks respectively, is the massive boson polariza- 
tion vector, / and /' are the Dirac spinors for the fermions and Vff is the unit matrix when 
//' are leptons and the CKM matrix when //' are quarks. The fcj and hi represent the 
structure of the vector-axial vector (V-A) coupling of the bosons, where for the case i = W 
we have kw = hw = l, i-e. purely left-handed coupling. Using the above coupling to 
fermions, it can be shown that the differential cross-section pp W/W — > ivX (figure |l]) 
for {W^ , W'^)^ is given by: 



d. GlM^ ^ I VP[^G;,(1 + z^) + 2AG~^,z] (2.2) 



drdi/dz 1927r 



qq' 



where z = cos 9 is defined as the scattering angle between the u-type quark and the outgoing 
neutrino (both being fermions) in the centre of mass (COM) frame, y is the rapidity of the 
intermediate boson, r = s/s is the ratio of the squares of the quark COM energy to the 
proton COM energy. S = S{s) and A = A{s) are functions of the quark COM energy. 

S= P^,{C^Cj)\CiCJy{k,kJ+h,hjf (2.3) 

i,j=W,W' i,j=W,W' 

A= Ai= Yl P^3{CiC,Y{CiC,f{k,hi + hjkif (2.4) 

i,j=W,W' ij=W,W' 



where 



, (s - Mf ){s - Mf) + TiTjMiMj 



where Mi, Fj are the mass and width of boson i respectively. The functions G^^, which 
appear in the differential cross-section are even or odd products of parton density functions 
for the relevant hadrons given by: 

^tg' " [lAixa, s)q'B{xb, s) ± qsixb, s)qA{xa, s)] (2.6) 

where qA,B{x,s) is the parton density function for a quark q carrying hadron momentum 
fraction a; in a collision in which the quark pair COM energy is s. The A, B indices represent 
the 'left' (travelling in the positive z-direction) or 'right' (travelling in the negative z- 
direction) hadrons respectively. This definition allows for easy modification of the pp — > 
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W/W' iuX cross-section to the pp — > W/W iuX, by changing the PDFs accordingly. 
Analogous expressions can also be written in the case of the (W~ ,W'~) by appropriately 
modifying the functions G^, and taking z —z. The width can be taken to be a free 
parameter in the reference model: the couplings of the W to other gauge bosons or the 
Higgs boson are model-dependent ^. We shall assume here for illustration that the fermionic 
decay width scales with the mass as Tw'^ff = (4rvi//3Mvi/)Miy/ (provided that M^y » 
Mt, the mass of the top quark) and that the tri-boson W'WZ vertex is suppressed by a 
small mixing angle and hence can be neglected in the analysis. 

2.2 W — W interference 

The narrow width approximation (NWA) is often used when discussing the production of 
new vector bosons. This approximation is usually claimed to be valid up to OiVw' /^w) 
corrections. But W-W interference effects can become important in certain regions even 
as the width T^/ 0, see for example [11], and as we also show below. Use of the NWA 
may thus lead to invalid conclusions, as pointed out in [19]. 

We expect to observe interference effects in the differential cross-section simply because 
the Drell-Yan process pp — > W/W' — > ivX can proceed either via a SM W or a W' in the 
reference model. The matrix element squared for the process may be decomposed in the 
following way: 

\M\^ = \Mw\^ + IMwl"^ + 2Re{M*wMw') (2.7) 

It is easy to see that the interference term depends on the functions S{s) and A{s) (eq. |2.3| 
and eq. |2.4| ). Here we discuss the function S{s) when studying interference effects, although 
the arguments for A{s) are equivalent. Just as with the squared matrix element, S{s) can 
be decomposed into pieces which are due to the W and W' individually and an interference 
piece: 

S = Sw,w + Sw',w' + Sw,w' + Sw',w = Sw,w + Sw',w' + 2Sint (2.8) 

where we have defined the interference term Sint = Sw,W' = Sw',w- Explicitly, this 
interference term may be written as: 



int 



{s - M^){s - M^,) + TwTw'MwMw 
[{s-M^r + Tl,M^][W^W'] 



{l + hwhw'f (2.9) 



where we have set all the couplings C^^^, = 1 and k\Y = kw/ = 1. It is evident that 
when hw = 1 (SM) and h^/ = — 1 then Sint = 0. Hence there is no interference for the 
case of the SM W and right-handed W' , and the square of the total matrix element for the 
process can be written as the sum of the squares of the individual matrix elements for the 
W and W': 

\M{hw' = = \Mwf + \Mw'\'^ (2.10) 

This is what we would expect just by noting that the W' decays to different final state 
particles than the W. However, when hw = 1 and h\Y' = 1, i.e. both left-handed, we have 
Sint 7^ 0. In fact, by examination of the expression for Sij (eq. ^), we can observe that Sint 



^An exception is the photon, for which the couphng is fixed by gauge invariance. 
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should be of the same order of magnitude as Sy/',w' ^-iid Sw,w- figure § shows the variation 
of the interference term for the case M^y/ = 1 TeV as well as Sw,w and Sw',W'- We can 
observe that Sint is negative in the intermediate mass squared region s G (M^, M^,). The 
total cross-section in this region is reduced in comparison to the sum of the individual W 
and W cross-sections. It is important to note that the interference term is non- vanishing 
and comparable in magnitude to the other terms in S{s) even as Vyi/i ^ 0, a clear indication 
of why NWA is not justified in the intermediate region. If Ty^i = the particles do not 
overlap directly with each other, however off-mass-shell effects still cause interference. We 
emphasize the fact that the interference is negative when a SM W interferes with a left- 
handed W . This leads to a reduction in the cross section with respect to the Standard 
Model expectation, a possibility seldomly considered in experimental searches. 



0.1 



0.00005 



-0.00005 



-0.1 




Figure 2: The interference term Sint in the case where hw' = 1, Mw' = 1 TeV, plotted against s, 
for different widths: Tw' = 1, 35, 100 GeV (green, blue, purple respectively). The terms Sw,w (red) 
and Sw' ,W' (black) are shown for comparison. It is evident that Sint is negative in the intermediate 
region (M^,M^,). It is also clear that the magnitude of the interference term is comparable to 
Sw,w and Sw',w'- As the width decreases the negative peak becomes narrower, but there always 
exists a portion of the curve which is independent of the width. 
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3. Extension to NLO 



Next, we extended the simulation to NLO using the MCSNLD and the POWHEG methods. 
The MCONLO method has previously been applied to the hadroproduction of gauge boson 
pairs [20,21], heavy quark-antiquark pairs [22] and single-top production [23]. In these 
appHcations, the Fortran Monte Carlo event generator HERWIG [24] was used to simulate 
the parton shower and hadronization. Within the framework of Herwig++, it has been 
applied to e+e~ annihilation to hadrons and Drell-Yan vector boson production, [25]. 

The method is based upon a careful expansion of the NLO results, in order to match 
certain features of the event generator used, in this case Herwig++. The shower approxi- 
mation to the NLO matrix element in Herwig++ must be subtracted from the exact NLO 
result in order to avoid double counting. This subtraction generates a number of negative 
weighted events which however, are few enough so that the number of events required for 
a smooth distribution is comparable to leading order simulations. 

The POWHEG method on the other hand generates the hardest emission of the parton 
shower to NLO accuracy first and, for angular ordered showers such as Herwig++, adds a 
truncated shower of soft and wide angled emissions between the hard scale and the scale 
of the hardest emission. The resulting partons are then showered subject to a veto so 
that no subsequent emissions have pt greater than the hardest emission. This method has 
been applied to Z pair production [26], heavy flavour production [27], Drell-Yan vector 
boson production [28,29], e^e~ annihilation into hadrons [30] and into top pairs and their 
decays [31], and NLO Higgs boson production via gluon fusion [32]. 

The advantage of this method over the MCSNLO method is the independence of the 
procedure on the event generator used and the generation of only positive weighted events. 

We now briefly discuss both methods and their application to W' boson production. 
Full details of the application of the MCONLO method to vector boson production can be 
found in Section 6 of [25]. Details of the application of the POWHEG method can be found 
in Chapter 4 of [33] where vector boson production is discussed in detail. The implemen- 
tation of a truncated shower of at most one emission in the framework of Herwig++ is also 
discussed. 



3.1 The MCONLO method 

The NLO cross-section for the production of W bosons can be written as a sum of two 
contributions, 

where Oqqi is the contribution from qcf annihilation and o'(q^q')g is the contribution from 
the Compton subprocesses. In the modified minimal subtraction (MSbar) factorization 
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scheme, these are: 



(To j dxidx2 



x[Dq{xi)Dg/{x2) + q'] 



Dq{Xq)Dqi{Xqi) 



5{1-x) + —Cf{-2 



27r 



l + x^ 
1-x 



Inx 



+ 4(1 +x- 



2 /ln(l - 



1 — X 



+ 



+ -TT^ ) <5(1 - X) 



ax\ax2 — — — r t^—Tr 



<i,<i' 

„2 



+ (x^ + (l+x^))ln 



2tt' 



- + 3x — -x 



(3.2) 



72 

where ao is the Born differential cross-section ^Q-f^y ^^^^ '5 the invariant mass and Y 
the rapidity of the vector boson. The xi,X2 are the NLO momentum fractions and Xq,Xq' 
are the Born momentum fractions with = XgXq'S, if S denotes the hadronic centrc-of- 
mass energy. Also, x = and Dq{xi) = xifq{xi) etc., with fq{xi) being the parton 

distribution function of parton q. 

Focusing on the qq' annihilation process for the moment, if we introduce the variable 



y = cos y. 



(3.3) 



where 6 is the scattering angle of the emitted parton in the partonic COM frame, we can 
re- write as an integral over x and y: 



^0 



dxdy 



x[Dq{xi)Dqt{x2) + q ^ q']l 

2 



Dq{Xq)Dq>{Xq>) 



"5, 



d{l -x) + ^Cf -2 



27r 



1 + x^ 



lux 



+ 4(1 



2 /ln(l - x] 



l-x 



+ (-8 + -7r2)<5(l-x; 



Mqq>{x,y) \+Mqq,{x,y) 



(3.4) 



where Mqq'{x,y) is the real emission matrix element. Since we have subtracted this con- 
tribution from the total cross-section, in the curly brackets we are left with the sum of the 
Born, virtual and QCD PDF correction contributions. Now we can define an infrared-safe 
observable O whose NLO expectation value is given by: 



{O''') = Y. dxdy 



Ow 
1 + 

1 — X 



■:[Dq{xi)Dq'{x2) + g ^ q] 1 
Dq{Xq)Dq,{Xq,) 2 

lnx + 4(l + x2)^^''^^~^) 



(5(1 - x) 



l-x 



+ 



+ (-8 + -7r2)<5(l-x; 



- Mqq, (x, y) \ + Ow'gMqq, (x, y) 



(3.5) 



where O^i and Oyi/ig are observables arising from hadronic final states generated from 
q + q ^ W and q + q ^ W + g starting configurations respectively. This however is not 
entirely correct because of double counting in the final states represented by Ojy/ which are 
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already included in the states arising from Owg- The solution to this is the subtraction 
of the parton shower contributions which we denote Mc^-,{x,y) from the regions in which 
the parton shower contributes (the jet region J) and integrate the full matrix element in 
the hard emission region D, left untouched by the shower. This gives for {O'^^'): 



{O 



qq 



E 



O 



w 



x[Dg{xi)Dg,{x2)+q^ q']l f as 



+ 4(1 + x 

+ Uw'n < M, 



Dq{Xq)Dq,{Xq. 

ln(l-x)\ / „ 2 



_l + x2 

-2- mx 



1 



8+^vr2)<5(l-x) 



+ E 



x[Dq{xi)Dqi{x2) + 9 ^ q'] 1 



D 



Ow 



+ 4(1 + 



+ Ow'gMqq,] 



Dq{Xq)Dq,{Xq>) 

ln(l-x)\ / . 2 



as 



-U(l-.) + ^C. 



_l + x2 

-2- mx 



X 



1 



+ ( -8 + |7r2 ) <5(1 - x) 



Mqq' 



(3.6) 



A similar functional {Oii^^')9) can be generated for the Compton subprocesses. Events can 
then be generated in the different regions of phase space according to their contributions 
to the above integrals. These events are then interfaced with Herwig++ and showered. Full 
details of the algorithm for event generation can be found in [25]. 

3.2 The PDWHEG method 

This method as described in [15] involves the generation of the hardest radiation from 
the parton shower according to the real emission matrix element and independently of the 
shower Monte Carlo generator used. If we introduce: 



R^^r = Mqqi + M, 



H.q,q')9 ' 



(3.7) 



where Mqq and M(^g^gi)g are real emission matrix elements for qq' annihilation and the 
Compton subprocesses respectively, we can write the cross section for the hardest gluon 
emission event as: 

da = Y^ Bld^^ [A«(0) + M^{pT)Rv,rd^r] ■ (3.8) 



The index q runs over all quarks and anti-quarks. The subscript v represents the Born 
variables, which in this case are the invariant mass Q and the rapidity Y of the boson, r 
represents the radiation variables x,y and d^y,d^r are the Born and real emission phase 
spaces respectively. 

A''(pt) is the modified Sudakov form factor for the hardest emission with transverse 



momentum pT, as indicated by the Heaviside function in the exponent of eq. 3.9 



A" (pa 



exp 



j d<^>rRv,r@{kTiv,r) - pt) 



(3.9) 



-8- 



where fcx is the transverse momentum of the hardest emission relative to the splitting axis 
and in this case is given by: 



4x 



Furthermore, 



(3.10) 



(3.11) 



Bv is the sum of the Born, Bv, virtual, Vv and real, R^^r terms, (with some counter-terms, 
Cv,r)- The Born variables are generated with distribution By, with the radiation variables 
of the first emission generated according to [A'^(O) + A'i{pT-)Ry^rd^r]- 

In the MSbar factorization scheme, the contribution to the order as cross-section for 



W production is given in eq. 3A and eq. 3^. The function B'^ in eq. 3.11 can then 
be written down as a sum of finite terms using the subtraction method. In this paper. 



we borrow the MC@NLO subtraction formula introduced in eq. 3.4 and write a function 
Bi{Q'^,Y) as: 



B^Q^^Y) = ^ / dxdydQ^dY 



dQ^dY 



x[Dq{xi)Dq>{x2) + q ^ q] 1 



1 + x^ 
-2- In a; +4(1 + x 



Dq[Xq)Dq:{Xq) 

2. ( ln(l — x) 



1 



1 



X 



+ 



3 



(5(1 -x) 



+ 



+ (x^ + (1 + x^))ln 



(1 -x)^ 



27r ^2 



+ 3x 



('l,<l')9 



+ M, 



{g,<i')g 



(3.12) 



where we have written the virtual and PDF corrections in terms of the real emission matrix 
elements and Mq are the subtracted parton shower approximation terms in the Herwig++ 
jet regions. Note that the above prescription does not imply that the POWHEG method 
depends on the shower MC used. We have simply used the shower approximation terms to 
defne a subtraction scheme for the definition of the NLO cross-section. 

The flavour of the event, the Born variables and Y, as well as radiation variables 



X and y are then generated according to the integrand in eq. 3.12. The radiation variables 
are ignored which amounts to integrating away these variables leaving the Born variables 
distributed according to B'^{Q^ ,Y). The radiation variables x,y are generated according 
to: 

A3{pT)R{x,y)dxdy . (3.13) 



Details of the algorithm used can be found in 
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4. Results 



We present a sample of distributions of variables obtained for ~ 10^ events using the 
Wpnlo event generator, both at leading and next-to-leading order, using the MCONLO (with 
Herwig++) and POWHEG methods. The general purpose event generator Herwig++, version 
2.2.1 [34], was used. The X-factor (where K = (Tnlo/^^lo) for the considered invariant 
mass range and for factorisation/renormalisation scales set to the default NLO scale = 
^ (where kj- and Q are the dilepton transverse momentum and invariant mass 
respectively) was found to be -ftT ~ 1.3, in all cases. The plots have been normalised to 
unity (apart from figure |lO|) to emphasise the differences in the shape of the distributions. 

For validation purposes, figure ^ presents a comparison of the W boson transverse 
momentum distribution, (assuming no W'^ between Tevatron data (taken from [35]) and 
the three possible methods: leading order, MC@NLO/Herwig++ and POWHEG. The plots include 
events in the invariant mass range (70 - 90) GeV. The MC@NLO/Herwig++ and POWHEG 
distributions are evidently in agreement with the data within the statistical Monte Carlo 
and experimental uncertanties. The leading order distribution is cut off at the W mass 
since this provides the only relevant scale in the shower, whereas the MC@NLO/Herwig++ 
and POWHEG distributions extend to higher transverse momentum. 

The subsequent figures in this section represent simulations made for the forthcoming 
CERN LHC running at 14 TeV proton-proton centre of mass energy, expected to run in 
the second half of 2009. figure shows the variation of the NLO cross section for a 1 TeV 
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Figure 3: Transverse momentum distribution at the Tevatron obtained for MC@NLO/Herwig++ in 
the DIS and MSbar factorisation schemes (PDFs: cteqSd and cteqSm [36] respectively), POWHEG 
MSbar (cteq5m) and LO (PDF: MRST2001LO [37]), in the mass range (70 - 90) GeV. 



left-handed W in the invariant mass range [400,5000] GeV with factorisation scale, [ip, for 
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Figure 4: The normalised variation witli scale of the cross section calculations at NLO (red 
crosses) and LO (black circles) are shown for a proton-proton collider at 14 TeV, AIw' = 1 TeV, 
Tw' = 36 GeV and left-handed chirality in the invariant mass range [400,5000] GeV. 



fixed renormalisation scale using the MSbar scheme. The LO variation with PDF scale is 
also shown in an equivalent range. The values have been normalised to the cross sections 



at the default scales /J-o = y k'^ + Q'^ (default NLO) and ^iq = Q (default LO). In the NLO 
case the renormalization scale was held fixed at My/i. The NLO cross section calculation 
appears to be slightly more stable than the LO calculation. The X-factor at /xq was found 
to be -fC = 1.288 and the LO cross section at /io = Q was found to be aio = (2.99±0.07)pb. 
Figures ^ and ^ show the transverse mass distributions at LO and NLO for a W at masses 
and widths of (1 TeV, 36 GeV) and (2 TeV, 72 GeV), for purely left-handed {hw = 1) 
and purely right-handed {hw' = — 1) couplings to fermions respectively. Figures |^ and ^ 
show the corresponding W/W' transverse momentum distributions. In this case the LO 
distribution cuts off at the W mass. The effect is less visible for higher W masses. 
Figure^ shows a comparison, for the different methods, of the W/W rapidity, longitudinal 
momentum and mass distributions for a right-handed W of mass 2 TeV and width 72 GeV 
at the LHC. 



Finally, figure 10 shows a comparison between the left- and right-handed transverse 
mass (Mr) distributions at NLO, using the MC@NLO/Herwig++ method. The importance of 
the interference between the SM W and the W' can be clearly observed: the differential 
cross-section in the region below Mt = M\^i in the purely left-handed case is reduced in 
comparison to the purely right-handed case. For transverse mass greater than the on-shell 
mass of the W , the interference term becomes positive for the left-handed case, although 
the effect is not significant. The SM contribution, in the absence of a W boson, is given 
for comparison in both figures. It can be observed that in the right-handed case the 
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contribution of the W' is simply additive to the SM contribution. 




Figure 5: Transverse mass distributions at the LHC obtained using the MC@NLO/Herwig++ and 
POWHEG methods (cteqSm/MSbar) and LO (MRST2001LO) for a purely left-handed W. The plots 
correspond to masses/widths equal to [1 TeV, 36 GeV] (left) and [2 TeV, 72 GeV] (right). The 
invariant mass range was taken to be (0.4 — 3.0) TeV for the 1 TeV case and (0.4 — 5.0) TeV for 
the 2 TeV case. The effect of the destructive interference can be observed in both cases. Note that 
the plots have been normalised to unity. 




Figure 6: Transverse mass distributions at the LHC obtained using the MC@NLO/Herwig++ and 
POWHEG methods (cteq5m/MSbar) and LO (MRST2001LO) for a purely right-handed W. The 
invariant mass range and W mass and widths are identical to those in the previous figure. 



5. Extraction of limits 

In Appendix ^ we provide a general method for discriminating between two models given 
a set of data. Here we apply this method to extract observation limits on the W mass and 
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Figure 7: Transverse momentum distributions at tlie LHC obtained using the MCSNLO/Herwig++ 
and POWHEG methods (cteqSm/MSbar) and LO (MRST2001LO) for a purely left-handed W. The 
invariant mass range and W mass and widths are identical to those in the previous figures. 




Figure 8: Transverse momentum distributions at the LHC obtained using the MCSNLO/Herwig++ 
and POWHEG methods (cteqSm/MSbar) and LO (MRST2001LO) for a purely right-handed W. The 
invariant mass range and W' mass and widths are identical to those in the previous figures. 



width at LO. A stand-alone program was written to calculate the quantity R at matrix 
element-level, given by eq. A. 5 : 

AT / irnw \ / a7 \ 



r^PiS) sr^ ( p{MtAT) \\ (Nt\ 



-{Nt-Ns) 



The 'true' underlying theory, called T, was taken to contain a W' at a pre-defined mass 
and theory S was taken to be the SM. Some comments are appropriate: 

• Although the total W width was being varied, the decay width to fermions was 
always taken to be Vw'^ff = {^Vw /'iMw)Mw' ~ 36 x Mw' {Mw' in GeV). 
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Figure 9: W/W rapidity (top left), longitudinal momentum (top right) and mass (bottom) dis- 
tributions at the LHC obtained using the MC(§NLD/Herwig++ and POWHEG methods (cteqSm/MSbar) 
and LO (MRST2001LO) for a purely right-handed W of mass 2 TeV and width 72 GeV. The 
invariant mass range and W mass and widths are identical to those in the previous figures. 



• In the experimental case the W mass would be unknown and maximum likelihood 
methods should be used to fit the parameters if significant deviation from the SM is 
found. 

• The R parameter can become very large if a small number of unlikely events occur, 
which favour one theory over the other. This is not realistic experimentally since 
unlikely events could arise from background or detector effects. To take into account 
these effects, one has to introduce nuisance parameters whose behaviour, at this 
level of analysis, has to be chosen arbitrarily. In the present analysis we avoid the 
introduction of such arbitrary parameters. 

• The exclusion curves were drawn for specific data distributions and fluctuations are 
expected. In other words, the plots given correspond to a specific experimental data 
set. 

• The ratio of the prior probability distributions, p{S)/p{T), was taken to be equal to 
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Figure 10: Transverse mass distributions at the LHC obtained using the MC@NLO/Herwig++ method 
(cteqSm/MSbar), POWHEG (cteqSm) for purely left- and right-handed W"s. The invariant mass range 
was taken to be (0.4 — 3.0) TeV. The plots correspond to masses/widths equal to (1 TeV, 36 GeV) 
and (2 TeV, 72 GeV). The significance of the destructive interference can be observed in the left- 
handed case; in the right-handed case the distribution is just the sum of the standard model W 
and right-handed W' contributions. Note that the plots are normalised to the NLO cross-section 
for each process. 



unity throughout this analysis: i.e. we assume both models are equally likely prior 
to the 'experiment'. 

• A rapidity cut on the leptons corresponding to Ucut = 2.5 for the LHC and ycut = 1-3 
for the Tevatron was applied to take into account the acceptance regions of the 
detectors. 

• The distributions p{Mt\S) and p{Mt\T) were calculated using the Monte Carlo event 
generator itself at higher statistics (~ 10^) than the required number of events to 
reduce the required computer time. The sum over i in eq. [A.5| was taken over the 
bins of these distributions and not individual events. 

The limits were drawn on a width-mass plane as log R = C (C is a constant) exclusion 
curves. The variable R can be interpreted as a probability ratio and an exclusion curve 
logR = C can be interpreted as the limit where the existence of a W is excluded with 
certainty 1 — . For example if C = 10, then the exclusion curve represents the ~ 
99.9996% confidence level. The LO exclusion curves can be seen, for different integrated 



luminosities at the LHC (14 TeV), in figure 11 for a right-handed W' and figure for a left 



handed W. The curves correspond to a single data sample at each {My^ri ^Ty^^i) point, and 
therefore there are large statistical fluctuations, particularly in the low-luminosity curves. 
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Figure 11: The detection reacti at ttie LHC for logi? = 10 (left) and log J? = 100 (right) at 
different integrated luminosities for the right-handed case. The colour scheme is; green, blue, black 
corresponding to the luminosities 1, 10, 100 fb^^. 
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Figure 12: The detection reach at the LHC for logi? = 10 (left) and logi? = 100 (right) at 
different integrated luminosities for the left-handed case. The colour scheme is identical to the 
previous figure. In the logi? = 10 and 100 fb^^ case all points below the contour have logi? < 10. 



A comparison between the curves for a left- and right-handed W is shown in figure 13. It 
can be observed that a left-handed W has a slightly higher detection reach, especially at 
higher widths. By examining figures |l^ and 12, we can deduce that the maximum detection 
reach at the LHC, for example using an integrated luminosity of 100 fb^^, for a W decaying 
primarily to fermions {Twi Tw'-^ff'), is ~ 4 TeV. We also show the expected limit at the 
Tevatron (1.96 TeV) in figure [l^ with an integrated luminosity of 2 fb~^, both at leading 
and next-to-leading (see below) orders. Note that the current experimental limit on the 
W' mass is 1 TeV, extracted from a sample of 1 fb~^ of data from the DO experiment [38]. 
When the W is only allowed to decay to fermions, i.e. has width Tiyi ~ 36 GeV, the 
predicted detection limit for logi? ~ 10 is Myyi 1.1 TeV. This is slightly better than 
the current Tevatron limit, but is expected to be reduced by experimental effects. Since 
the available centre-of-mass energy at the Tevatron is 1.96 TeV, we expect the saturation 
of the detection reach to come at about Myy ~ 1 TeV without interference and slightly 
higher in the left-handed case when interference effects are included. The Tevatron NLO 
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Figure 13: The detection reacli at tlie LHC for logi? = 10 (left) and logi? 100 (right) at 
different integrated luminosities for the left- and right-handed cases. The colour scheme is for 1, 
10, 100 fb~^ is: left-handed: green, blue, black and right-handed: orange, light blue, pink. 



case does not exhibit any substantial difference from the LO case. 
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Figure 14: The detection reach at the Tevatron for logi? 10 (left) and logi? — 100 (right) at 2 
fb"^, for the left- and right-handed cases, at LO and NLO. The colour scheme is for right-handed 
and left-handed correspondingly, at LO: green, blue and NLO: light blue, pink. 



We have performed an equivalent analysis using the NLO method POWHEG at matrix 
element level to improve computational time. Working at matrix element level with the 
POWHEG method is justified since the transverse mass distribution is not significantly altered 
after shower and hadronisation and no difficulties arise due to negative-weighted events, as 
would be the case in the MCONLO case. The comments given at the beginning of the section 
for the LO analysis also apply to the NLO analysis. The results are shown in figures |l5| 
and |l^ in comparison to the LO results. In the right-handed chirality scenario, NLO implies 
a lower detection reach than indicated at LO. The situation is more complicated in the 
left-handed case where the NLO case implies a slightly higher reach for larger widths. 

To investigate the dependence of the NLO results on the factorisation scale fip we 
have reproduced the log R contours for the right-handed W LHC case with an integrated 
luminosity of 10 fb~^ at different values oi while keeping the normalisation scale fixed, 
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Figure 15: The detection reacti at ttie LHC for logi? = 10 (left) and logJi" = 100 (rigiit) at 
different integrated luminosities for the right-handed case compared at LO and NLO. The colour 
scheme for 1, 10, 100 fb^^ is: LO: green, blue, black and NLO: orange, light blue, pink. 



logR = 10 



logR = 100 




2000 2500 3000 3500 



4500 5000 5500 




M (GeV) 



2000 2500 3000 3500 4000 4500 5000 5500 



M (GeV) 



Figure 16: The detection reach at the LHC for logi? = 10 (left) and logi? = 100 (right) at 
different integrated luminosities for the left-handed case compared at LO and NLO. The colour 
scheme is identical to the previous figure. 



using the MSbar scheme. The results are shown in figure 17. The curves show that 
the factorisation scale does not affect the detection reach substantially, for example only 
shifting the log R = 10 contour at a width of Tw' ~ 200 GeV from Mw' ~ 3500 GeV to 
Mw' ~ 3750 GeV going from f^p = 0.5/io to fip = AfiQ. 



6. Conclusions 

We have presented a Monte Carlo implementation of the Drell-Yan production of new 
charged heavy vector bosons. We have considered the interference effects with the Stan- 
dard Model W boson, allowing arbitrary chiral couplings to the leptons and quarks. 
Moreover, the implementation is correct up to next-to-leading order in QCD, via the 
MC@NLO/Herwig++ and POWHEG methods. We have presented a sample of results at both 
leading and next-to-leading orders. As expected, the LO and NLO boson transverse mo- 
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Figure 17: The NLO detection reach at the LHC for logi? = 10 (left) and logi? = 100 (right) 
for an integrated luminosity of 10 fb~^ at different factorisation scale fip- The colour scheme for 
Hf = 0.5/io, Mo, 2/j,o and 4/io is: green, blue, pink and light blue. 

mentum distributions were found to differ significantly, the NLO extending to higher pT- 
The dilepton transverse mass, invariant mass, rapidity and z-momentum distributions were 
found not to be significantly altered by the NLO treatment. The total cross section was 
found to increase in the NLO case by a factor of ~ 1.3 in the region of interest. 

Subsequently we applied a theoretical discrimination method to the W' reference model 
to obtain mass-width observation curves for left- and right-handed chiralities of the W 
both at LO and NLO (POWHEG). The NLO curves were shown not to vary significantly 
with factorisation scale. The event generator used throughout this analysis, Wpnlo, is fully 
customisable and publicly available [10]. 
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A. Model discrimination 

The search for new physics is essentially a task of discriminating between two models: 
one with new physics, the other without. The actual task of finding any new physics 
though is laborious: one has to understand the detector well enough and has to be able to 
obtain enough statistical significance to say with certainty that something new has been 
discovered. Here we adopt a rather theoretical approach: we describe a purely statistical 
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method for discriminating between models [39,40]. This will essentially yield an upper 
bound on the detection reach of a heavy particle: detector effects and backgrounds will 
result in a reduced detection limit. It is useful, however, to be aware of the theoretical 
possibilities for discovery. 

A.l Likelihood ratios of probability density functions 

Consider N data points, of a mass variable measurement, {rrii}. Based on these data 
points, a theoretical model T is i? times more likely than another theoretical model S, if, 

'^ = 1^ (A.l) 

p{S\{mi}) 

where p{X\{mi}) is the probability of model X being true given the data set {m,}. We 
may use Bayes' Theorem to rewrite R as 

p{{m}\S)p{S) 

where p{T) and p{S) are the probabilities that S and T are true respectively, usually called 
prior probabilities since they represent previous knowledge on the theories. We assume 
that these quantities are equal: there is no strong evidence for either model. We may 



simplify eq. A. 2 further: 



PjS) ^ IltoP{m\T) ^ ^ p{m,\T) 



where we have assumed that events in the data set {rrij} are independent, eq. A.3 is a 
discrete version of the Kullback-Leibler distance [41], a useful quantity for comparing the 
relative likelihood of two theories according to a data sample. However, it is important to 
note that the distributions p{mi\T) and p{mi\S) are normalized to unity. This means that 
any difference in the number of events predicted by the two theories will not be taken into 
account. This will obviously underestimate the significance of a difference in number of 
events, for example a substantial excess of events in an invariant mass peak. We describe 
a method which takes this factor into account in the next section. 

A. 2 Poisson likelihood ratios 

In this modification to the method described in the previous section, we simply multiply 



the variable R defined in eq. AT by a ratio of Poisson distributions for the total number 
of events: 



N 



p{T\{mi}) ( Nt\ ^_(^Mt-Ns) 



^= o ^ e-^''^-'''> (A.4) 

piS\{mi}) \NsJ 

where Nx = ctx-L is the expectation value of the number of events according to theory 
X, given by the product of the cross-section, ax, and the integrated luminosity, L. This 
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expression can be manipulated in a similar manner to eq. A.3 to obtain: 



N 



^1=0 



p{m,\S)J \Ns' 



(Nt-Ns) 



(A.5) 



For convenience we may define the 'shape' and 'Poisson' factors respectively: 

'p{mi\T)' 



N 



Rs = exp ^ log 



i=0 



p{mi\S) 



Rp 



N 



Ns ' 



(A.6) 



This method takes into account the difference in the total number of events expected 
according to each theory at the given integrated luminosity. This is accomplished by re- 
weighting the 'shape' factor Rs by a factor Rp which gives the ratio of probabilities to 
obtain the observed number of events. 

A. 2.1 Application to a toy model 




Figure 18: The differential cross-sections ^ according to two 'toy' theories T and S arc shown. 
Theory T possesses a Gaussian 'bump', at m = 5 whereas S is just a falling distribution. 1/m. m 
is in arbitrary mass units and a in equivalent inverse area squared units. 

Before applying the method to the full W model, it is instructive to present its appli- 
cation to a simple model involving two analytical 'toy' distributions. Events for the two 
distributions have been generated by the general Monte Carlo event generation method. 
The 'differential cross-sections' for the two 'theories' T and S with respect to a variable m 
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(in arbitrary units) are given by (defined in the range [0.1, 10]): 



dm m 
das 1 



+ 0.3e-^™-°^ (A.7) 

(A.8) 



dm m 

Theory T has a Gaussian peak at m = 5 on top of a background falhng as ~ 1/m and 



theory S falls as ~ 1/m. The situation is shown in figure 18. This is qualitatively similar 
to the SM tail (theory S) and the SM plus a heavy particle (theory T). The 'cross-sections' 
in the range m = [0.1, 10] were found to be gt = 5.14 and as = 4.60, in arbitrary area 
units. Assuming an integrated 'luminosity' of L = 30 (equivalent arbitrary inverse area 
units), we have an expected number of events Nt = 154 and Ns = 138. We assume that 
theory T is the correct underlying theory: events are actually distributed according to it. 
The result for the variable R if theory T was 'true' was found to be i? = 62. This implies 
that theory T is 62 times more likely than theory S given the specific data set. If, however, 
the underlying theory is chosen to be S, then R = 0.23. Note that in the case that theory 
T is 'true', it is easier to exclude theory S than to exclude theory T in the case that theory 
S is 'true'. 
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